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Abstract 

We consider a measurement of exclusive production of scalar Xc(0 ++ ) meson in the proton- 
proton collisions at LHC and RHIC and in the proton-antiproton collisions at the Tevatron via 
XcO — > 7r + 7r~ decay. The corresponding amplitude for exclusive double-diffractive XcO meson pro- 
duction was obtained within the /cj-factorization approach including virtualities of active gluons 
and the corresponding cross section is calculated with unintegrated gluon distribution functions 
(UGDFs) known from the literature. The four-body pp —> ppir + 7r~ reaction constitutes an irre- 
ducible background to the exclusive XcO meson production. We calculate several differential distri- 
butions for pp(p) —> pp(p)XcO process including absorptive corrections. The influence of kinematical 
cuts on the signal-to-background ratio is investigated. Corresponding experimental consequences 
are discussed. 
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I. INTRODUCTION 



The mechanism of exclusive production of mesons at high energies became recently a 
very active field of research (see e.g. Ref. [l[ and references therein). Central exclusive 
production processes represent a very promising and novel way to study QCD in hadron- 
hadron collisions. Recently, there is a growing interest in understanding exclusive three-body 
reactions pp — > ppM at high energies, where the meson (resonance) M is produced in the 
central rapidity region. In particular, these reactions provide a valuable tool to investigate in 
detail the properties of resonance states. Many of the resonances decay into tttt and/or KK 
channels. The representative examples are: M = a, p°, /o(980), 0, /2(1275), /o(1500) and, as 
we shall emphasis here, Xco- Various decay channels can be studied. It is clear that these 
resonances are seen (or will be seen) "on" the background of a tttt or KK continuum 1 . The 
two-pion background to exclusive production of / (1500) meson was already discussed in 
Ref. Q. The recent works concentrated on the production of Xc mesons (see e.g. Refs. (HMol 
and references therein) where the QCD mechanism is similar to the exclusive production 
of the Higgs boson. Furthermore, the Xc(o,2) states are expected to annihilate via two- 
gluon processes into light mesons and may, therefore, allow the study of glueball production 
dynamics (lb! ]. 

Recently, the CDF Collaboration has measured the cross section for exclusive production 



of Xc mesons in proton-antiproton collisions at the Tevatron [11], by selecting events with 
large rapidity gaps separating the centrally produced state from the dissociation products of 
the incoming protons. In this experiment Xc mesons are identified via decay to the J/ip ±7 
with J ftp — > n + p~ channel. The experimental invariant mass resolution was not sufficient 
to distinguish between scalar, axial and tensor Xc- While the branching fractions to this 
channel for axial and tensor mesons are large [l2[ (B = (34.4 ±1.5)% and B = (19.5 ±0.8)%, 
resp ectively) the branching fraction for the scalar meson is very small B = (1.16 ± 0.08)% 
[l2j . On the other hand, the cross section for exclusive Xco production obtained within the 
k t - factorization is much bigger than that for Xci and Xc2- As a consequence, all Xc mesons 
give similar contributions [7J to the J/ift + 7 decay channel. Clearly, the measurement via 
decay to the J/ip + 7 channel cannot provide cross section for different Xc- 

Could other decay channels be used? The Xco meson decays into several two-body chan- 
nels (e.g. tttt, K + K~, pp) or four-body hadronic modes (e.g. tt + tt~tt + tt~ , tt + tt~K + K~). 
The branching ratios are shown in Table HI In this paper we analyze a possibility to measure 
Xco via its decay to tt + tt~ channel. The advantage of this channel is that the tt + tt~ contin- 
uum has been studied recently [l3[ and is relatively well known. In addition the axial xa 
does not decay to the tttt channel and the branching ratio for the Xc2 decay into two pions is 
smaller. A much smaller cross section for Xc2 production means, in practice, that only Xco 
will contribute to the signal. 

In the present paper, we wish to calculate differential distributions for the exclusive 
production Xco meson with a few UGDFs taken from the literature relevant for small gluon 
virtualities (transverse momenta). We shall use matrix element for the off-shell gluons 
as obtained in Ref. 0. The expected non-resonant background can be modeled using a 
"non-perturbative" framework, mediated by pomeron-pomeron fusion with an intermediate 
off-shell pion exchanged between the final-state particle pairs. Thus, we consider pp{p) — > 



1 In general, the resonance and continuum contributions may interfere. This may produce even a dip. A 
good example is the /o(980) production (see e.g. Ref. @, Q). 
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TABLE I: Branching fractions for the Xcj two- and four-body hadronic decays, taken from Ref. 
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Channel 


B(Xco) 


B(Xci) 


B(Xc2) 


7T + 7T~ 

K+K- 
PP 

ir+ir-K+K- 


(0.56 ± 0.03) x 10" 2 
(0.610 ± 0.035) x 10~ 2 
(2.28 ±0.13) x 10" 4 
(2.27 ±0.19) x 10~ 2 
(1.80 ±0.15) x 10" 2 


(0.73 ± 0.04) x 10~ 4 
(0.76 ± 0.26) x 10~ 2 
(0.45 ±0.10) x 10" 2 


(0.16 ±0.01) x 10~ 2 
(0.109 ± 0.008) x 10~ 2 
(0.72 ± 0.04) x 10~ 4 
(1.11 ±0.11) x 10~ 2 
(0.92 ±0.11) x 10~ 2 



pp(p)7T + 7T~ reaction as a genuine four-body process with exact kinematics which can be 
easily used when kinematical cuts have to be improved. 

Exclusive charmonium decays have been a subject of interest at the e + e~ colliders as they 
are an excellent laboratory for studying quark-gluon dynamics at relatively low energies. 
Thus, a measurement of many exclusive hadronic Xc decays if possible is very valuable. 
Although these Xc states are not directly produced in e + e~ collisions, they are copiously 
produced in the radiative decays ip{2S) — > 7X C > each of which has a branching ratio of around 
9% [12|| . The CLEO Collaboration has studied exclusive ^ c (o,i,2) decays to four-hadron final 



7T 1 7T 



7r°7r°, 



K + K~ 



7T 7r°, pp7l°7T°, 



states involving two charged and two neutral mesons [14 
K + K~r]TT° and K ± 7T Zf K°7i°. The BESIII Collaboration has studied two-body Xc(o,2) decays 
into 7r°7r° and r/r/ 15] 2 and four-body Xc(o,i,2) decays into 7r 7r 7r 7r [l6[ final states, where 
Xcj signals appear in radiative photon energy spectrum. Recently (Ref. [l7[) a study of Xcj 
hadronic decays and measurements of Xcj A(1520)A(1520) decaying into ppK + K~ has 
been presented. In the present paper we start the discussion of the possibility to measure 
the different decay channels in proton-(anti)proton collisions in order to determine the cross 
section for exclusive production of the P-wave charmonia. Here, continuum backgrounds can 
be larger than in the e + e~ collisions and this requires a separate discussion of the feasibility. 
We will discuss this issue in the present paper. 

The paper is organized as follows. In section [TT] we give general expressions for the 
amplitudes of the exclusive Xco meson production and 7r + 7r~ pairs production including 
discussion of absorptive corrections. Section [TTT1 contains the presentation of the main results 
and a discussion of the uncertainties related to the approximations made. Finally, some 
concluding remarks and outlook are given in Section IIVI 



II. SIGNAL AND BACKGROUND AMPLITUDES 



A. Diffractive QCD amplitude for exclusive Xco production 



The QCD mechanism for the diffractive production of heavy central system has been 
proposed by Khoze, Martin and Ryskin (KMR) and developed in collaboration with Kaidalov 
and Stirling for Higgs production (see e.g. Refs. [l8-2l]|). In the framework of this approach 



2 The Xci decay into these final states are not considered as they are forbidden by the spin-parity conser- 
vation. 
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>XcO 



FIG. 1: The QCD mechanism of exclusive diffractive production of XcO meson including the ab- 
sorptive correction. 

the amplitude of the exclusive pp — > ppXco process is described by the diagram shown in 
Fig. [H where the hard subprocess g*g* — > Xco is initiated by the fusion of two off-shell gluons 
and the soft part represented in terms of the off-diagonal unintegrated gluon distributions 
(UGDFs) . The formalism used to calculate the exclusive Xco meson production is explained 
in detail elsewhere 0] and so we will only review relevant aspects here. 

The "full" amplitude for the exclusive process pp — > ppXco can be written as 



M f ^ l X PPXc (s,y,-p ljt ,-p 2tt ) 



= MliZppxX 3 ' y> -pi,» -P2, t ) 
+ M r ;xt Xc (s iyi -p htl -p 2 j). 



We can write the "bare" amplitude |5| as 



off, 



r 2Nf^j d 



(2.1) 



(2.2) 



where the objects fg{f g2 are skewed (or off-diagonal) unintegrated gluon distributions of 
both nucleons. t\ )2 are the momentum transfers along each nucleon line 3 , qij, q 2 ,t, %,t, %i,2, 
x[ ~ x' 2 <C X12 are the transverse momenta and the longitudinal momentum fractions for 
active and screening gluons, respectively. UGDFs are nondiagonal both in the x and q 2 
space. The usual off-diagonal gluon distributions are nondiagonal only in x. In the limit 
^1,2 - > x'i >2 , <?o,t — 1i/2t anc ^ ^1,2 ~~ they become the usual UGDFs. 

The vertex factor V(q lt ,q 2t ) describes the coupling of two virtual gluons to Xco meson 
is obtained in heavy quark approximation and can be written as 



V(q ltt , <?2,i) = k nlo 



8tg 2 n'{0) 3M 2 q ht q 2 j-2ql t ql t -(q ltt q 2tt )(ql t + ql i 



M y/nMN c 



(M 2 + q\, 



a 2 ) 2 

H2,t) 



(2.3) 



where M is the Xco mass, g 2 = 47ra s (M 2 ) and the strong coupling constant is calculated in 
the leading order and extended to the nonperturbative region according to Shirkov-Solovtsov 



In the following for brevity we shall use notation t% 2 which means t\ or t^.- 
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analytical model 22|. The value of the P-wave radial wave function at the origin is taken 
to be [23|| 7?.' (0) = V0.075 GeV 5//2 and the radiative corrections factor in the vertex K NLO 
is well-known [24| . K NLO ~ 1.68. 

The rescattering correction shown in Fig. [1] by the extra blob can be written in the form 



M 



rescatt 
pp^rppXc 



's,y,-Pi, t ,-P2,t) = J 



2(2tt) 



M 



bare 

pp->ppXc 



[s,y, fci, t ,fe 2 , 



t , 



(2.4) 



where t = — p 1 1 — k t and fc 2 ,t = — p 2t — fc t with momentum transfer k t . The amplitude 
for elastic proton-proton scattering at an appropriate energy is conveniently parameterized 
as 



A pp (s, k 2 ) 



Aq(s) exp{-Bk 2 /2) . 



(2.5) 



From the optical theorem we have lmA (s,t = 0) = scr tot (s) (the real part is small in the 
high energy limit). B is the effective slope of the elastic differential cross section 



B(s) = Bi + 2a'M 



(2.6) 



and is adjusted to the existing experimental data for the elastic NN scattering. The 
Donnachie-Landshoff parametrization |25] of the total pp or pp cross sections can be used 
to calculate the rescattering amplitude. We take s = 1 GeV 2 , a'p = 0.25 GeV" 2 and Bf. 



B% N = 9 GeV" 2 . 



The KMR UGDFs, unintegrated over q t , are calculated from the conventional (integrated) 
distributions g(x, qf) and the so-called Sudakov form factor \/T g (qf, /i 2 ) as follows: 



f« MR (x } x> J Ql^;t) = R, 



d 



d In g 2 



F(t). 



(2.7) 



9?= 



The Sudakov factor suppresses real emissions from the active gluon during the evolution, so 
that the rapidity gaps survive. The factor R g approximately accounts for the single InQ 2 
skewed effect [26j. In the calculations presented here we take R g = 1.3, and value of the 
hard scale is /1 2 = M 2 /4. The choice of the scale is somewhat arbitrary, and consequences 
of this choice were discussed in Ref. jBf. 

In the original KMR approach the following prescription for the effective transverse mo- 
mentum is taken: 



Q\ t = min (gg t , q 2 l t ) , Q\ t = min (g 2 t , g 2 J . 



{21 



Other prescriptions are also possible jB], @|. In the KMR approach only one effective gluon 
transverse momentum is taken explicitly in their skewed UGDFs compares to two indepen- 
dent transverse momenta in our case (see Eq. (12.21) ). Please note also that the skewed KMR 
UGDFs does not explicitly depend on x' (assuming x' < a; C 1). 

In Ref. [27j a procedure was presented which allows to calculate the generalized (or 
skewed) parton distributions of the proton, H(x, £; g 2 , /i 2 ), unintegrated over the partonic 
transverse momenta, from the conventional parton distributions, q(x,fi 2 ) and g(x,/j, 2 ), for 
small values of the skewedness parameter ( 2 < 1 and any x. The momentum fractions 
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carried by the emitted and absorbed partons are defined as = x ± ^£ with support 
— 1 < x < 1. The result is a simple approximate phenomenological form for the distribution: 



ding 2 ' 2?r ft • -~^ + q 2 



dxgjx, q 2 ) 2 iV c a 3 / \x + /i 2 
^ oiTj + % ) -T— In — + 1.2- 



+ 5 — 1 (a^tvai^, g t 2 ) + a;(i va i(x, q 2 )) . (2.9) 

Z7T 



x=5 



In evaluating / g 's we have used the GRV NLO [28] and GJR NLO |29| collinear gluon 
distributions, which allow to use rather low values of gluon transverse momenta Q\ = g 2 4 , 
Qiti Q2,t — 0-5 GeV 2 . The collinear distributions such as CTEQ and MRST are defined 
for higher factorization scales (Q 2 > 1 GeV 2 ), and therefore are less useful in applications 
discussed here. 

The i-dependence of the unintegrated gluon distribution f g 's is not well known and is 
isolated in the effective form factors of the QCD Pomeron-proton vertex, which are param- 
eterized in the forward scattering limit in the exponential form as 

F(t) = exp(6t/2) (2.10) 

with the t-slope parameter 6 = 4 GeV -2 . Then the integral in Eq. (12.41) can be evaluated as 

11 



v, -Pi, t , -P2, t ) = 47rg( ff ° + 26 ) exp 2(B + 2b) ) M %Z™°( 8 > y > ~ Pl ^ ~ P2 ^ ■ 

(2.11) 

The cross section for the three-body reaction is calculated as 

a = j -\M\ 2 (2n) 6 (p B + ft - Pl - ft - ft ) w {27r)32E3 (2.12) 

by choosing convenient kinematical variables. 

B. Diffractive amplitude for 7T+7T - continuum 

The dominant mechanism of the exclusive production of 7r + 7r~ pairs at high energies is 
sketched in Fig. |2j In calculations of the amplitude related to double diffractive mechanism 
for the pp — > pp7r + 7i~ reaction we follow the general rules of Pumplin and Henyey j3l| used 
recently in Ref. [l3j 4 . 

The "full" amplitude for the exclusive process pp — > ppn + 7i^ (with four-momenta p a + 
Pb — > Pi + T>2 + P3 + Pi) can be written formally as 



M£!*(s,y»,y*,p„,ih.,P„) = M bar ' 



pp-*-pp7T7rl S ) U3i Vii Pl,t) P2,t) Pm) ~ M pp^ppmrK 8 ; > 1/3) 2/4, Pi,*) P-2p Pm) 

pp-^ppirni 3 'i 2/3, 2/4, Pl,£, P2,ti Pm) , 



(2.13) 



4 For early rough estimates see Ref. [32[. There have also been a variety of experimental results, in particular 
from the CERN ISR (for recent reviews see [lj). 
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- 7T (Pi) 




FIG. 2: The double-diffractive mechanism of exclusive production of 7r + 7r pairs including the 
absorptive corrections. 

where the auxiliary quantity p m = p 3 t — p it . 
The "bare" amplitude can be written as 



M 



bare 

pp—tppnir 



M 13 {t 1 ,s 13 )F 7V {t a ) 
+ M u (t 1 ,s 14 )F 7r (t b 



mi 



1 



h - ml 



-F n {t a )M 24 {t 2 ,S 24 ) 

F w (t b )M 23 (t 2 ,s 23 ) , 



(2.14) 



where denotes "interaction" between nucleon i — 1 (forward nucleon) or % = 2 (backward 
nucleon) and one of the two pions k = 3 (7r + ), k = 4 The energy dependence of the 

amplitudes of the nN subsystems was parameterized in the Regge form with pomeron and 
reggeon exchanges. The details of the matrix element are explained in Ref. [l3|. The 
strength parameters and values of the pomeron and reggeon trajectories are taken from the 
Donnachie-Landshoff analysis [25] of the total cross section for ttN scattering. The slope 
parameters of the elastic nN scattering can be written as shown in Eq. ( 12.61) and are Bp 
= 5.5 GeV~ 2 , = 4 GeV -2 [l3[, for pomeron and reggeon exchanges, respectively. 

The Donnachie-Landshoff parametrization (25| can be used only for the nN subsystem 
energy > 2 — 3 GeV (see Ref. [IH). Bellow = 2 GeV the resonance states in 
re N subsystems are present. In principle, their contribution could and should be included 
explicitly 5 . In order to exclude resonance regions we shall "correct" the parametrization 
Eq. (12.141) ) multiplying it by a purely phenomenological smooth cut-off correction factor 



eX p (IMo) 



f nN (w \ — [ i 

JconA W ik) - — 7 

1 + exp I 



w ik -w \ 

I / 



(2.15) 



The parameter W 



2 GeV gives the position of the cut, and the parameter a = 0.2 GeV 

^ont^Wik) ~ 1 an d close to kmematical 

~ m n + m^) ss 0. 

The extra form factors F(t a ) and F(tb) "correct" for off-shellness of the intermediate 
pions in the middle of the diagrams shown in Fig. [2j In the following they are parameterized 



describes how sharp the cut-off is. For large energies fcmt(Wik 
threshold ftZtiW* 



The higher the center-of-mass energy the smaller the relative resonance contribution. 
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as 

F«(taj>) = exp ( %F^*) , (2-16) 

i.e. normalized to unity on the pion- mass-shell. In general, the parameter A a ff t E is not 
known precisely but, in principle, could be fitted to the (normalized) experimental data. 
From our general experience in hadronic physics we expect A 2 Q ^ E ~ 1.2 — 2 GeV 2 . How to 
extract A //,.e will be discussed in the result section. 

The absorptive corrections to the "bare" amplitude (Eq. (12.141) ) can be written as: 



Vf''''"'"' 1 / d2k t A Pp( S > k t) Kji,„n , , ,,-)-■ 

Jvl pp-+ppmr — 1 / 2(27r) 2 e J 1 PP^PPtttt\P(:.! Pi. r Ph. I P-2.il- ( - ■ I M 



reseat ■ I u "** ^ppl" 5 )"^ \ Abare 



where p* = p a - fc t , p b * = p 6 + fc t . 

The amplitude described above (see Eq. 12. 13f) is used to calculate the corresponding cross 
section including limitations of the four-body phase-space. The cross section for the two-pion 
continuum is obtained by integration over the four-body phase space: 



d 3 pi d 3 p 2 d 3 p 3 d 3 p4 
(2tt) 3 2E 1 (2ir) 3 2E 2 (2vr) 3 2£ 3 (27r) 3 2£ 4 ' 

(2.18) 



In order to calculate the total cross section one has to calculate the eight-dimensional integral 
numerically. The details how to conveniently reduce the number of kinematical integration 
variables are given elsewhere [HI]. 



III. RESULTS 

Let us start from presenting various differential cross sections. In practical integrations of 
the exclusive Xco meson we choose the transverse momenta of outgoing nucleons (j>i,t,P2,t), 
the meson rapidity (y) and the relative azimuthal angle between outgoing nucleons (^12)- 

In Fig. [3] we show distributions of the central exclusive Xco production cross section at 
y/s = 14 TeV without (dashed line) and with (solid line) absorptive corrections as described 
in subsection III Al These calculations were done with GJR NLO j29[ collinear gluon dis- 
tribution, to generate the KMR UGDFs (see Eq. (12. 7j) ). which allows to use low values of 
the gluon transverse momenta Q\ > Q 2 cut = 0.5 GeV 2 . The bigger the value of the cut-off 
parameter, the smaller the cross section (see Ref. |5j). In the calculations we take the value 
of the hard scale to be fi 2 = M 2 /4. The smaller fi 2 , the bigger the cross section jBJ. The 
absorption effects lead to a damping of the cross section. In most distributions the shape 
is almost unchanged. Exception is the distribution in proton transverse momentum where 
the absorption effects lead to a damping of the cross section at small proton p t and an en- 
hancement of the cross section at large proton p t . In relative azimuthal angle distribution 
we observe a "diffractive dip" structure in the region of <p%2 ~ tt/2. Transverse momen- 
tum distribution of Xco shows a small minimum at p t ~ 2.5 GeV. The main reason of its 
appearance is the functional dependence of matrix elements on its arguments j5[. 

In Fig. H] we compare distributions of the central exclusive Xco production cross section 
at yfs = 14 TeV calculated for two collinear gluon distributions: GRV94 NLO (upper lines) 
and GJR08 NLO (bottom lines). We show results here only for distributions with absorptive 
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FIG. 3: Differential cross sections for the pp — > ppXcO reaction at y/s = 14 TeV without (dashed 
line) and with (solid line) absorption effects. These calculations were done with the GJR08 NLO 



29( UGDFs. 



corrections calculated with the KMR off-diagonal UGDFs given by Eq. (12. 7p (solid lines) and 
with off-diagonal UGDFs given in the phenomenological form given by Eq. (I2.9p . The peaks 
at large rapidities appear only when we use formula (12. 7p . In this region one of off-diagonal 
UGDFs changes a sign. This shows limitations in applying formula (12. 7\i . 

Let us turn now to an estimation of the two-pion background. In Fig. [5] we show the two- 
pion invariant mass distribution at the center-of-mass energy of the CERN ISR yfs = 62 
GeV (this is the highest energy at which experimental data exist). Here we have used 
a simple model described in subsection III Bl The experimental cuts on the rapidity of 
both pions and on longitudinal momentum fractions (Feynman-x, xp = 2p\\/y/s) of both 
outgoing protons are included when comparing our results with existing experimental data 
from 133] (see Refs. (34], |35[ for early studies). Absorption effects were included in this 



calculation. The results depend on the value of the nonperturbative, a priori unknown 
parameter of the form factor responsible for off-shell effects (see Eq. (I2.16P ). We show 
results with the cut-off parameters ^-off,E = 1-8, 2.0 GeV 2 as represented by the dotted, 
dashed and solid line, respectively. The experimental data show some peaks above our flat 
model continuum. They correspond to the well known 7r + 7r~ resonances: a, p°, ^2(1275) 
which are not included explicitly in our calculation. In the present analysis we are interested 
mostly what happens above > 3 GeV, above the resonance region where there are no 
experimental data points. Our model with A. ff,E parameter fitted to the data provides an 
educated extrapolation to the unmeasured region. 

In Fig. [6] we show differential distributions in pion rapidity for the pp — > ppn + n~ reaction 
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FIG. 4: Differential cross sections for the pp — > ppXcO reaction at y/s = 14 TeV with absorption 
effects. The results with KMR off-diagonal UGDFs given by Eq. (|2.7h (solid lines) and with off- 
diagonal UGDFs given by Eq. (12. 9p (dashed lines) are shown. These calculations were done for 
two UGDFs: GRV94 NLO (upper lines) and GJR08 NLO (bottom lines). 




FIG. 5: Differential cross section da/dM n7T for the pp 



ppir+n 



reaction at y/s = 62 GeV with 



experimental cuts relevant for the CERN ISR experimental data from Ref. 33J. Results with the 
cut-off parameters ^-offE = 1-6,1.8,2.0 GeV 2 are shown by the dotted, dashed and solid line, 
respectively. The absorption effects were included in the calculations. 



at y/s = 0.5, 1.96, 14 TeV without (upper lines) and with (bottom lines) absorption effects. 
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The integrated cross section slowly rises with incident energy. The camel-like shape of the 
distributions is due to the interference of components in the amplitude: pomeron-pomeron 
component dominates at midrapidities of pions and pomeron-reggeon (reggeon-pomeron) 
peaks at backward (forward) pion rapidities, respectively (see Ref. [13(). The reader is 
asked to notice that the energy dependence of the cross section at y n ~ is reversed by the 
absorption effects which are stronger at higher energies. 




FIG. 6: Differential cross section da/dy n for the pp — > ppir + Tr~ reaction at yfs = 0.5, 1.96, 14 TeV 
with A-ojfE = 2.0 GeV 2 . The results without (upper lines) and with (bottom lines) absorption 
effects are shown. 

Now we wish to compare differential distributions of pions from the Xco decay with those 
for the continuum pions. In the first step we calculate the two-dimensional distribution 
da(y,pt)/dydp t , where y is rapidity and p t is the transverse momentum of Xc- The decay of 
Xco —> 7r + 7r _ is included then in a simple Monte Carlo program assuming isotropic decay of 
the scalar Xco meson in its rest frame. The kinematical variables of pions are transformed 
to the overall center-of-mass frame where extra cuts are imposed. Including the simple cuts 
we construct several differential distributions in different kinematical variables. In Fig. [7] 
we show distributions in pion transverse momenta. The pions from the decay are placed at 
slightly larger transverse momenta. This can be therefore used to get rid of the bulk of the 
continuum by imposing an extra cut on the pion transverse momenta. 

In Fig. |H] we show two-pion invariant mass distribution for the double-diffractive nir 
continuum and the contribution from the decay of the Xco meson (see the peak at M n7T ~ 3.4 
GeV). In these figures the resonant Xco distribution was parameterized in the non-relativistic 
Breit-Wigner form: 

^(XcO -> 7T + 7T~) a pp ^ ppXc0 2M n7r -j—p, ... =rz , (3.1) 



dM^ ^ ^pp^ppxco 7T (M 2 W — M 2 ) 2 + (M n7T T 



with parameters according to PDG In calculation of the Xco distribution we use GRV94 
NLO and GJR08 NLO collinear gluon distributions. The cross sections for the Xco production 
and for the background include absorption effects. While upper row shows the cross section 
integrated over the full phase space at different energies, the lower row shows results including 
the relevant pion pseudorapidities restrictions — 1 < T] w +,r] w - < 1 (RHIC and Tevatron) and 
-2.5 < r] n+ ,7] n - < 2.5 (LHC). 
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FIG. 7: Differential cross section do jdpt^ at yfs = 0.5, 1.96, 14 TeV with cuts on the pion pseu- 
dorapidities. The diffractive background was obtained with the cut-off parameters A.X £ = l-6> 
2.0 GeV 2 (lower and upper dashed lines, respectively). Results for the pions from the decay of the 
XcO meson including the ir + ir~ branching ratio, for GRV94 NLO (upper lines) and GJR08 NLO 
(bottom lines) UGDFs, are shown. The absorption effects were included in the calculations. 
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FIG. 8: The 7r + 7r~ invariant mass distribution at \fs = 0.5,1.96,14 TeV integrated over the full 
phase space (upper row) and with the detector limitations in pion pseudorapidities (lower row). 
The dashed lines present the -kit continuum with the cut-off parameters A^j E = 1.6, 2.0 GeV 2 
(lower and upper dashed lines, respectively). The XcO contribution is calculated with GRV94 NLO 
(dotted line) and GJR08 NLO (solid line) collinear gluon distributions. The absorption effects for 
the XcO meson and for the background were included in the calculations. 

The question now is whether the situation can be improved by imposing extra cuts. In 
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Fig. [9] we show results with additional cuts on both pion transverse momenta \pt,w\ > 1-5 
GeV. Now the signal-to-background ratio is somewhat improved especially at the Tevatron 
and LHC energies. Shown are only purely theoretical predictions. In reality the situation is, 
however, somewhat worse as both protons and, in particular, pion pairs are measured with a 
certain precision which leads to an extra smearing in M n7T . While the smearing is negligible 
for the background, it leads to a modification of the Breit-Wigner peak for the Xco meson 
6 . The results with more modern GJR UGDF are smaller by about factor of 2-3 than those 
for somewhat older GRV UGDF. 
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FIG. 9: The ir + ir~ invariant mass distribution at y/s = 0.5, 1.96, 14 TeV with the relevant restric- 
tions in the pion pseudorapidities and pion transverse momenta. The dashed lines present the 
7T7T continuum with the cut-off parameters E = 1.6, 2.0 GeV 2 (lower and upper dashed lines, 
respectively). In calculating the Xco contribution we use GRV94 NLO (dotted line) and GJR08 
NLO (solid line) collinear gluon distributions. The absorption effects for the Xco meson and for 
the background were included. Clear Xco signal with relatively small background for the Tevatron 
and LHC energies can be observed when imposing extra cuts on p t ^ . 



In Table [TT] we have collected the numerical values of the cross sections (see a pp ^ PPXc0 in 
Eq. (13.11) ) for exclusive Xco production for some selected UGDFs at different energies. 



TABLE II: Integrated cross sections in nb (with absorption corrections) for exclusive Xco production 
at different energies with GRV94 NLO [3] and GJR08 NLO [2^] collinear gluon distributions. In 
this calculations we have taken the relevant limitations in the pion pseudorapidities 1%) < 1 at 
the RHIC and Tevatron, j^l < 2.5 at the LHC and lower cut on both pion transverse momenta 
|pt, ff | > 1.5 GeV. 





full phase space 


with cuts on r] n 


with cuts 


on ij v and p t;7r 


(TeV) 


GRV 


GJR 


GRV 


GJR 


GRV 


GJR 


0.5 


63.6 


34.5 


14.2 


7.9 


7.3 




4.0 


1.96 


309.0 


127.1 


51.1 


21.1 


25.8 




10.6 


14 


1152.3 


329.7 


554.1 


159.1 


183.7 




52.2 



6 An additional experimental resolution not included here can be taken into account by an extra convolution 
of the Breit-Wigner shape with an additional Gaussian function. 
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The main experimental task is to measure the distributions in the Xco rapidity and trans- 
verse momentum. Can one recover such distributions based on the measured ones in spite 
of the severe cuts on pion kinematical variables? In Fig. [10] we show the two-dimensional 
ratio of the cross sections for the Xco meson in its rapidity and transverse momentum: 

Ration, ft) = wi ] hcuts n ' V Pt . (3.2) 

The numerator includes limitations on r] n and p tj7r . These distributions provide a fairly 
precise evaluation of the expected acceptances when experimental cuts are imposed. The 
experimental data could be corrected by our two-dimensional acceptance function to recover 
the distributions of interest. 





FIG. 10: Ratio of the two-dimensional cross sections in (y,pt) for the pp —> ppXcO reaction with 
the relevant limitations on the pion pseudorapidities and a few lower cuts on the pion transverse 
momenta pt^- These calculations were done with GJR08 NLO 29] UGDFs. The upper row is for 
the STAR detector (\rj n \ < 1) and the lower row for the ATLAS or CMS detectors (1^1 < 2.5). 



IV. CONCLUSIONS 



It was realized recently that the measurement of exclusive production of Xc via decay 
in the J/ip + 7 channel cannot give production cross sections for different species of Xc- 
In this decay channel the contributions of Xc mesons with different spins are similar and 
experimental resolution is not sufficient to distinguish them. 
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In the present paper we have analyzed a possibility to measure the exclusive production 
of Xco meson in the proton-(anti)proton collisions at the LHC, Tevatron and RHIC via 
Xco — > 7r + 7r~ decay channel. Since the cross section for exclusive Xco production is much 
larger than that for Xci an d Xc2 an d the branching fraction to the mr channel for Xco is larger 
than that for Xc2 (Xci does not decay into two pions) the two-pion channel should provide 
an useful information about the Xco exclusive production. 

We have performed detailed studies of several differential distributions and demonstrated 
how to impose extra cuts in order to improve the signal-to-background ratio. The two- 
pion background was calculated in a simple model with parameters adjusted to low energy 
data. We have shown that relevant measurements at Tevatron and LHC are possible. At 
RHIC the signal-to-background ratio is much worse but measurements should be possible 
as well. Imposing cuts distorts the original distributions for Xco m rapidity and transverse 
momentum. We have demonstrated how to recover the original distributions and presented 
the correction functions for some typical experimental situations. 

In the present paper we have concentrated on 7r + 7r - channel only. This is mainly due to 
the relatively good control of the background. Measurements of other decay channels, e.g. 
K + K~, are possible as well and will be discussed elsewhere. 

Acknowledgments 

This study was partially supported by the Polish grant of MNiSW N N202 249235. 



M.G. Albrow, T.D. Coughlin and J.R. Forshaw, Prog. Part. Nucl. Phys. 65 (2010) 149. 

P. Lebiedowicz, A. Szczurek and R. Kamiriski, Phys. Lett. B680 (2009) 459. 

D. Aide et al. [GAMS Collaboration], Phys. Lett. 397 (1997) 350. 

A. Szczurek and P. Lebiedowicz, Nucl. Phys. A826 (2009) 101. 

R.S. Pasechnik, A. Szczurek and O.V. Teryaev, Phys. Rev. D78 (2008) 014007. 

R.S. Pasechnik, A. Szczurek and O.V. Teryaev, Phys. Lett. B680 (2009) 62. 

R.S. Pasechnik, A. Szczurek and O.V. Teryaev, Phys. Rev. D81 (2010) 034024. 

L.A. Harland-Lang, V.A. Khoze, M.G. Ryskin and W.J. Stirling, Eur. Phys. J. C65 (2010) 

433. 

L.A. Harland-Lang, V.A. Khoze, M.G. Ryskin and W.J. Stirling, Eur. Phys. J. C71 (2011) 
1545. 

C. Amsler and F.E. Close, Phys. Rev. D53 (1996) 295. 

T. Aaltonen et al. [CDF Collaboration], Phys. Rev. Lett. 102 (2009) 242001. 

K. Nakamura et al. [Particle Data Group], J. Phys. G37 (2010) 075021. 

P. Lebiedowicz and A. Szczurek, Phys. Rev. D81 (2010) 036003. 

Q. He et al. [CLEO Collaboration], Phys. Rev. D78 (2008) 092004. 

M. Ablikim et al. [BESIII Collaboration], Phys. Rev. D81 (2010) 052005. 

M. Ablikim et al. [BESIII Collaboration], Phys. Rev. D83 (2011) 012006. 

M. Ablikim et al. [BESIII Collaboration], arXiv:hep-ex/1103.2661. 

V.A. Khoze, A.D. Martin and M.G. Ryskin, Phys. Lett. B401 (1997) 330; 

V.A. Khoze, A.D. Martin and M.G. Ryskin, Eur. Phys. J. C23 (2002) 311. 

A.B. Kaidalov, V.A. Khoze, A.D. Martin and M.G. Ryskin, Eur. Phys. J. C33 (2004) 261. 

A.B. Kaidalov, V.A. Khoze, A.D. Martin and M.G. Ryskin, Eur. Phys. J. C31 (2003) 387. 



15 



[21] V.A. Khoze, A.D. Martin, M.G. Ryskin and W.J. Stirling, Eur. Phys. J. C35 (2004) 211. 
[22] D.V. Shirkov and I.L. Solovtsov, Phys. Rev. Lett. 79 (1997) 1209. 
[23] E.J. Eichten and C. Quigg, Phys. Rev. D52 (1995) 1726. 

[24] R. Barbieri, M. Caffo, R. Gatto and E. Remiddi, Nucl. Phys. B192 (1981) 61; W. Kwong, 
PB. Mackenzie, R. Rosenfeld and J.L. Rosner, Phys. Rev. D37 (1988) 3210; M.L. Mangano 
and A. Patrelli, Phys. Lett. B352 (1995) 445. 

[25] A. Donnachie and PV. Landshoff, Phys. Lett. B296 (1992) 227. 

[26] A.G. Shuvaev, K.J. Golec-Biernat, A.D. Martin and M.G. Ryskin, Phys. Rev. D60 (1999) 
014015. 

[27] A.D. Martin and M.G. Ryskin, Phys. Rev. D64 (2001) 094017. 

[28] M. Gliick, E. Reya and A. Vogt, Z. Phys. C67 (1995) 433. 

[29] M. Gliick, D. Jimenez-Delgado and E. Reya, Eur. Phys. J. C53 (2008) 355; 

M. Gliick, D. Jimenez-Delgado, E. Reya and C. Schuck, Phys. Lett. B664 (2008) 133. 
[30] V.A. Khoze, A.D. Martin and M.G. Ryskin, Eur. Phys. J. C24 (2002) 581. 
[31] J. Pumplin and F. S. Henyey, Nucl. Phys. B117 (1976) 377. 

[32] Y.I. Azimov, V.A. Khoze, E.M. Levin and M.G. Ryskin, Sov. J. Nucl. Phys. 21 (1975) 215. 

[33] A. Breakstone et al. (ABCDHW Collaboration), Z. Phys. C48 (1990) 569. 

[34] R. Waldi, K.R. Schubert and K. Winter, Z. Phys. C18 (1983) 301. 

[35] A. Breakstone et al. (ABCDHW Collaboration), Z. Phys. C42 (1989) 387. 



16 



